Light-induced multistability and Preedericksz transition in nematic liquid crystals 
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We study light transmission tlirougli a tiomeotropically oriented nematic liquid crystal cell and 
solve self-consistently a nonlinear equation for the nematic director coupled to Maxwell's equations. 
We demonstrate that above a certain threshold of the input light intensity, the liquid- crystal cell 
changes abruptly its optical properties due to the light-induced Freedericksz transition, demonstrat- 
ing multistable hysteresis-like dependencies in the transmission. We suggest that these properties 
can be employed for tunable all-optical switching photonic devices based on liquid crystals. 
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I. INTRODUCTION 

Liquid crystals (LCs) play an important role in the 
modern technologies being used for numerous applica- 
tions in electronic imaging, display manufacturing, and 
optoelectronics ili, i2j ■ A large variety of electro-optical ef- 
fects that may occur in LCs can be employed for a design 
of photonic devices. For example, the property of LCs to 
change its orientational structure and the refractive in- 
dex in the presence of a static electric field suggests one 
of the most attractive and practical schemes for tuning 
the photonic bandgap devices 0, Ql • Nonlinear optical 
properties of LCs and multistability of light transmission 
are of a great interest for the future applications of LCs 
in photonics pj- 

Light polarized perpendicular to the LC director 
changes its orientation provided the light intensity 
exceeds some threshold value Q- This eflFect is 
widely known as the light-induced Freedericksz transition 
(LIFT), and its theory was developed more than two 
decades ago in a number of the pioneering papers [3,lal3- 
In particular, Zeldovich et al. J\ demonstrated that the 
light-induced Freedericksz transition can generally be 
treated as the second-order orientational transition, but 
in some types of LCs hysteresis-like dependencies and 
two thresholds can be observed, for the increasing and 
decreasing intensity of the input light. The results ob- 
tained later by Ong [l^l confirmed that for the MBBA 
nematics the Freedericksz transition is of the second or- 
der and there is no hysteresis behavior, whereas for the 
PAA nematics the Freedericksz transition is of the first 
order and the hysteresis-like behavior with two distinct 
thresholds should be observed. Although these conclu- 
sions have been confirmed to some extent in later experi- 
ments [l^ , the theory developed earlier was based on the 
geometrical optics and by its nature is approximate. The 
similar approximation was used later jllj for taking into 
account a backward wave in a LC film placed in a Fabry- 
Perot resonator, where it was shown that the threshold 
of the Freedericksz transition depends periodically on the 
LC cell thickness. 

Nonlinear optical properties of a nematic LC film in 



a Fabry-Perot interferometer was studied by Khoo et 
al. Il2| , who considered the propagation of light polarized 
under an acute angle to the LC director and observed 
experimentally bistability in the output light intensity 
caused by giant nonlinearity of the LC film. Cheung et 
al. 13] observed experimentally the effects of multista- 
bility in a similar system, including oscillations of the 
output light intensity. 

However, in spite of numerous theoretical studies and 
experimental observations, a self-consistent theory of the 
light-induced Freedericksz transition based on a system- 
atic analysis of the coupled equations for the nematic 
director and electromagnetic field is still missing. There- 
fore, the purpose of this paper is twofold. First, we 
consider a general problem of the light transmission 
through a homeotropically-oriented nematic LC and an- 
alyze the specific conditions for the multistability and 
light-induced Freedericksz transition, for possible appli- 
cations in all-optical switching photonic devices. Sec- 
ond, for the first time to our knowledge, we consider this 
problem self-consistently and solve numerically a coupled 
system of the stationary equations for the director and 
Maxwell's equations. We present our results for two kinds 
of nematic liquid crystal, para-azoxyanisole (PAA) and 
Np-methoxybenzylidene-np-butylaniline (MBBA) , which 
show quite dissimilar behavior of the nematic director 
at the Freedericksz transition in the previous theoreti- 
cal studies [l3, and also discuss light transmission and 
bistability thresholds as functions of the cell thickness. 

The paper is organized as follows. Sections II and III 
present our basic equations and outline our numerical 
approach. Section IV summarizes our results for two 
kinds of nematic liquid crystal and discusses in detail 
both bistability and hysteresis-type behavior of the light 
transmission. Section V concludes the paper. 



II. BASIC EQUATIONS 

We consider a nematic LC cell confined between two 
planes (z = and z — L) with the director initially ori- 
ented along the z axis (see Fig.P). The LC cell interacts 
with a normally incident monochromatic electromagnetic 



wave described by the electric field E(r, t), 

E(r, t) = i [E(r)e-*'^* + E*(r)e''^*] . (1) 



Je = --—eikEiEl , Eik = e^Sik + SaniUk- 

OTT 



E, 




FIG. 1: (colour online) Schematic representation of the prob- 
lem. A LC cell is placed between two walls (z — 0,z = L), 
the vector n describes the molecules orientation in the cell. 



Here /d is the LC elastic energy density, fs is a contri- 
bution to the free energy density from the light field, n 
is the ncmatic director, Ku are the elastic constants, Eik 
is the LC dielectric permittivity tensor, Eq = e|| — e^ > 
describes anisotropy of the LC dielectric susceptibility, 
where ey and e± are the main components of the tensor 
Eik parallel and perpendicular to the director, respec- 
tively. 

We assume that outside the LC cell the electric field 
is directed along the x axis (see Fig.QJ, which can cause 
the director reorientation in the xz plane inside the LC 
cell. When the incident beam is broad, we can describe 
it as a plane wave, so that all functions inside the LC 
cell will depend only on the z-coordinate. Therefore, we 
can seek the spatial distribution of the nematic director 
in the form 



To derive the basic equations, we write the free energy 
of the LC cell in the presence of the electromagnetic wave 
as follows |3 



F^jfU 



fE)dV, 
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[n X V X n]^, 



n(r) = ex sin 0(z) + e^ cos 4>{z), (3) 

where (f) is the angle between the director and the z axis 
(see Fig. ^ , ex and e^ are the unit vectors in the Carte- 
sian coordinate frame. 

After minimizing the free energy |(2Jl with respect to 
the director angle (j), we obtain the stationary equation 
for the LC director orientation in the presence of the light 
field 
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where we take into account that, as follows from 
Maxwell's equations, the electric vector of the light 
field inside the LC cell has the longitudinal component 
E.iz) = -{E,x/e,,)Exiz). 

From Maxwell's equations, we obtain the scalar equa- 
tion for the z-component of the electric field. 



d^Ex 
dz^ 



E I E 



±fc|| 



e_L + £a COS^ 



tEx = 0, 



(5) 



where k = 2ttX/c, and A is the wavelength of the incident 
light. The time- averaged ^-component of the Poynting 
vector, Sz = (c/8'k)ExH* remains unchanged inside the 
LC cell 13, [l3 , and it can be used for characterizing dif- 
ferent regimes of the nonlinear transmission. 



III. NUMERICAL APPROACH 

We solve the system of coupled nonlinear equations 
Q and |(SJ) in a self-consistent manner together with the 
proper boundary conditions. For the director, we assume 
the strong anchoring at the cell boundaries, i.e. 



0(0) = 0(L) = 0, 



(6) 



whereas for the electric field we consider the standard 
scattering conditions 



Ex{0) — Ein + Erci, Ex{L) — Eo 



(7) 



Here ii^in, i?rcf , and ii^out are the amplitudes of the inci- 
dent, reflected, and outgoing waves, respectively. In all 



equations above we consider the magnetic susceptibility 
/i = +1, and the refractive index outside the cell rig — 1, 
also taking into account that Hy — {l/ik){dEx/dz). 

The boundary conditions |2J imply that we consider 
two counter-propagating waves on the left side of the LC 
cell, incoming and reflecting, whereas only an outgoing 
wave appears on the right side. Therefore, in order to 
solve this nonlinear problem, first we fix the amplitude 
of the outgoing wave Eout ■ It allows us to find the unique 
values of the incident Ein and reflected E^cf waves. 

Equation for the director Q is similar to a general- 
type equation for a nonlinear pendulum with the fixed 
boundary conditions ©. This means that we should 
look for its periodic solutions with the period 2L. In fact, 
there exist many periodic solutions of Eq. 10}. First of 
all, a trivial solution ^(z) — corresponds to the undis- 
turbed orientation distribution of the director and the 
absolute minimum of the free energy Q- The Freeder- 
icksz transition occurs when this trivial solution becomes 
unstable for larger values of the input light intensity, and 
the director angle (l){z) becomes nonzero. We find this 
solution numerically by using the well-known shooting 
method pl4| . By fixing the amplitude of the outgoing 
wave i?out and taking (j){L) = at the right boundary, 
we find the values of the derivative {d(j) / dz) z=l such that 
after integration we obtain a vanishing value of the direc- 
tor angle at the left boundary, (f){Q) = 0. By analyzing 
the nonlinear equation Q in a two-dimensional phase 
space, we can show that the corresponding solution lies 
just below the separatrix curve, and it has no node be- 
tween the points z = and z ^ L. This observation 
allows us to reduce significantly the parameter region for 
the required values of the derivative (dip/ dz)z=L- From 
the obtained set of solutions we chose the solution that 
corresponds to the absolute minimum of the free energy 

We also take into account the fact that a finite energy 
barrier can appear between the minima of the free energy 
which correspond to the trivial and nontrivial solutions 
for the director orientation angle (j){z). When the light 
intensity decreases adiabatically, the director does not 
return to its initial undisturbed position at the threshold 
value of the " up" Freedericksz transition, but it remains 
in a disturbed state which corresponds to a local mini- 
mum of the free energy; thus, the transition to the state 
4>{z) = takes place only when this energy barrier disap- 
pears. This leads to a hysteresis-like dependence of the 
director and the different threshold values for the "up" 
and "down" transitions in the director orientation. 



IV. RESULTS AND DISCUSSIONS 

We solve the nonlinear transmission problem for two 
kinds of nematic liquid crystals, para-azoxyanisole (PAA) 
and Np-methoxybenzylidene-np-butylaniline (MBBA), 
which possess different signs of the parameter B = 
(1 - 9e||/(4ej^) - [K^^ - Kn) / K^^i) / A, which appears in 
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FIG. 2: (colour online) (a,b) Spatial distributions of the field 
amplitude \Ex\ in the cell of MBBA, before (dashed) and after 
(solid) the light-induced Freedericksz transition, L = A/no, 
A = 6328 A, no = 1.544. (b) Spatial distributions of the 
director deviation angle (fiiz) in the cell of MBBA after the 
light-induced Freedericksz transition, for L = A/no (solid), 
L = 100/im (dashed), are shown together with the function 
(fiQ sin{n z / L) at (po = 1.483 (dash-dotted). 



the geometrical optics approximation |2l ll(l | . According 
to their approach , the sign of this parameter B deter- 
mines the order of the Freedericksz transition. For PAA 
B < 00 and the Freedericksz transition should be of the 
first order, while for MBBA B > and there should be 
the second order transition. 

We take the following physical parameters [l3|: (a) 
for PAA, Kii = 9.26 ■ IQ-'^dyn, K33 = 18 • IQ-'^dyn, 
no ^ 1.595, He == 1.995, at A = 4800 A, and (b) for 
MBBA, Kii = 6.95 • IQ-'^dyn, K33 = 8.99 • I0'''dyn, 
no = 1.544, Tie = 1.758, at A = 6328^; and consider two 
values for the cell thickness, L — X/uq and L — lOO/im. 

Spatial distributions of the electric field amplitude 
£'2,(2;) I in the LC cell before and after the light-induced 
Freedericksz transition occurs is presented in Fig. |2ta) 
for the parameters of MBBA and the cell thickness 
L = X/uq. For the other value of the LC cell thickness 
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FIG. 3: (colour online) Maximum deformation angle (^max vs. Sz in PAA for: (a) L — A/no, (b) L — lOO/im, at A = 4800^1 
and no = 1.595, and in MBBA for: (c) L = A/no, (d) L = lOO^m, at A = 6328 A, and no = 1.544. 



(L — lOO^m), the spatial distribution of the electric field 
is similar, but the number of the oscillations of the elec- 
tric field \Ea:\ inside the LC cell increases due to a larger 
value of L/X. For PAA, a very similar distribution of 
the electric field is found. Thus, we reveal an essentially 
inhomogeneous spatial distribution of the electric field 
inside the LC cell, and the functions |i?2;(z)| are different 
before and after the Freedericksz transition. 

Spatial distribution of the director orientation angle 
(j){z) inside the LC cell after the Freedericksz transition 
is shown in Fig. Efb) for the parameters of MBBA, for 
L = A/no and L = lOO/um, respectively. On the same 
plot, we show the function (po sm{Trz/L) at (po = 1.483 for 
comparison. We notice that the position of the maximum 
of the director deviation angle can shift from the point 
z = L/2, as a consequence of an asymmetric distribution 
of the field |i?a;(z)| inside the LC cell. Spatial distribution 
of the director angle (j>{z) in the PAA cell has the same 
character as that shown in Fig. Efb) for MBBA. 

In Fig. 121 we present our numerical results for a change 
of the maximum deformation angle 0max of the director 



as a function of the power density inside LC Sz for in- 
creasing and decreasing light intensity, for both PAA and 
MBBA and two values of the cell thickness, L = A/no 
and L = lOOjim. For both kinds of LC, we observe a 
hysteresis-like dependence of the angle ^max and two dif- 
ferent thresholds of the light-induced director reorienta- 
tion: S'z, for the increasing intensity, and Sz , for the 
decreasing intensity. In both cases, these two thresh- 
olds correspond to the first-order transition. The results 
are similar for two values of the LC cell thickness, see 
Figs.|3|[a-d). Thus, our results suggest that at the light- 
induced Freedericksz transition the cells of both kinds 
of LCs, PAA and MBBA, reveal hysteresis-like behavior 
with the respect to Sz- 

Dependencies of the amplitude of the outgoing wave 
l^-outl on the amplitude of the incident wave |i?in| are 
shown in Figs. 0Ja-d) for the parameters of both PAA 
and MBBA. Depending on the LC cell thickness L, the 
cell transmission is characterized by either hysteresis or 
multistability with respect to the incident wave ampli- 
tude. In the case of small thickness of the LC cell 
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FIG. 4: (colour online) Multistable transmission of the LC cell, shown as the outgoing wave |i?out| vs. the incident wave \Ei„ 
for PAA: (a) L = A/no and (b) L = lOOfim, and for MBBA: (c) L = A/no and (d) L — lOOfim, respectively. 



(L — A/rio) only the hysteresis-like transmission is ob- 
served; it is caused by the hysteresis behavior of the direc- 
tor reorientation between "up" and "down" thresholds, 
as presented in Figs.0Ja,c). However, for larger thickness 
(L = lOOfxm) we observe the transmission multistability, 
above the "up" threshold for increasing light intensity, 
and above the "down" threshold for decreasing light in- 
tensity [see Figs.0Jb,d)]. Multistability in our system is 
similar to that of a nonlinear resonator, and is it deter- 
mined by the resonator properties of a finite thickness of 
the LC cell. 

The thresholds of the director reorientation for in- 
creasing and decreasing light intensities are shown in 
Figs.[SIa,b), for PAA and MBBA, respectively, as func- 
tions of the normalized thickness of the LC cell. Simi- 
lar to the results of the geometrical optics approxima- 
tion 13, [l3) the threshold values are proportional to 
(1/L)^, but they increase approximately in two times 
due to an essentially inhomogeneous spatial distribution 
of the electric field inside the LC cell. A similar increase 
of the threshold value for an inhomogeneous distribution 



of the electric field in the LC cell was also mentioned 
by Lednei et al. [l5|. In addition, for the "up" thresh- 
old we observe an additional periodic dependence of the 
threshold value on the cell thickness L, which is typical 
for resonators and is caused by an interference of two 
counter-propagating waves in the LC cell. This result 
agrees with the results obtained for LC in a Fabry-Perot 
resonator i^. The "up" threshold is determined by a 
competition between the electric field forces and elastic 
forces of the liquid crystal, and thus the interference dis- 
tribution of the electric field in the LC cell is important. 
However, the "down" threshold is defined by the condi- 
tion of the disappearance of a barrier between the local 
and absolute minima of the LC free energy [7| . We sup- 
pose that difference of these mechanisms leads to the dif- 
ferent type of L-dependencies for the "up" and "down" 
thresholds. 

We should mention that our results differ qualitatively 
from the results of earlier studies [3,ll3i where for MBBA 
both hysteresis and bistability were not predicted. In the 
simplest case of one traveling wave 0, [l3 j the conserva- 
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FIG. 5: (colour online) Thresholds of the director reorienta- 
tion for increasing (solid) and decreasing (dotted) light inten- 
sities vs. the cell thickness L: (a) PAA, (b) MBBA. 



tion of the value of Sz during the Freedericksz transition 
leads to the conservation of the electromagnetic field am- 
plitudes at the boundaries of the LC cell. However, in 
the general case there always exists a reflected wave, so 
that we have S^ = {c/8tt)E^{0)H*{0) = (c/87r)(^i„ + 

^rof)(^in-i;rof) = (c/Stt) (£;?„- £;2^f) = 5in-^rof, whcrC 

5'in, ^ref are the power densities of the incident and re- 
flected waves, respectively. In such a case, the conserva- 



tion of 5*2 does not require the conservation of S'in and 
S'ref , so that the amplitudes of the electromagnetic fields 
at z = can change at the Freedericksz transition, as is 
seen in Fig. |21 Thus, the problem solved in this paper 
and that in Refs. [3,[i3 corresponds to different bound- 
ary conditions. Therefore, we suggest that experimen- 
tally observed the second order Freedericksz transition 
for MBBA liquid crystal [a llg is caused by the weak 
reflection from the boundaries of LC cell. In that sit- 
uation, the single wave approximation can be used and 
results obtained in Refs. 0, U^ become valid. 
V. CONCLUSIONS 

We have analyzed the light transmission through 
homeotropically-oriented cell of a nematic liquid crystal, 
and studied multistability and light-induced Freedericksz 
transition. We have solved numerically the coupled sta- 
tionary equations for the nematic director and electric 
field of the propagating electromagnetic wave, for two 
kinds of liquid crystals (PAA and MBBA). We have found 
that the liquid crystals of both kinds possess multistabil- 
ity and hysteresis behavior in the transmission charac- 
terized by two thresholds of the director reorientation, 
so that for the increasing and decreasing light intensities 
the Freedericksz transition is of the first order. 

We have demonstrated that the resonator effects of the 
liquid-crystal cell associated with the light reflection from 
two boundaries are significant, and they arc responsible, 
in particular, for the observed periodic dependence of 
the threshold values and multistability of the transmit- 
ted light as a function of the cell thickness. We expect 
that these features will become important for the study 
of perio dic p hotonic structures with holes filled in liquid 
crystals [I3 where multiple reflection effects and nonlin- 
ear light-induced Freedericksz transition should be taken 
into account for developing tunable all-optical switching 
devices based on the structure with liquid crystals. 
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